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SUMMARY 


A  computerized  algorithm  to  generate  cross-sectional  dimensions  and  fiber 
orientations  for  composite  airframe  structures  is  described,  and  its  application 
in  a  wing  structural  synthesis  is  established.  The  algorithm  unifies  computa¬ 
tions  of  aeroelastic  loads,  stresses,  and  deflections,  as  well  as  optimal  struc¬ 
tural  sizing  and  fiber  orientations  in  an  open-ended  system  of  integrated  computer 
programs.  A  finite-element  analysis  and  a  mathematical-optimization  technique 
are  the  main  components  of  the  procedure.  Design  constraints  include  stress, 
strain,  local  buckling,  and  minimum  gage.  The  algorithm  is  applied  to  evaluate 
a  series  of  basic  design  alternatives  such  as  types  of  construction,  types  of 
material,  and  manufacturing  restrictions  for  the  low-aspect-ratio  wing  of  a  large 
transport  airplane.  Included  also  is  a  review  of  the  computational  efficiency 
of  the  method . 


INTRODUCTION 

Extensive  use  of  computer-aided  design  in  synthesis  of  airframes  for  a  low- 
aspect-ratio  supersonic  cruise  airplane  is  reported  in  references  1,  2,  and  3. 

The  studies  in  reference  1  focused  on  the  arrow-wing  configuration.  The  objec¬ 
tive  was  to  size  the  wing  structure  for  minimum  mass,  subject  to  static  strength 
requirements.  The  basic  wing  was  assumed  to  be  constructed  of  ribs,  spars,  and 
sandwich  cover  panels.  The  variants  considered  in  the  sizing  studies  were  an 
all-metal  (titanium)  wing,  a  wing  of  hybrid  construction  in  which  the  cover  panels 
were  made  of  a  graphite-polyimide  composite,  and,  finally,  an  all-composite  wing. 

Because  the  airplane  (ref.  1)  was  large  and  very  flexible,  it  was  important 
to  include  the  effect  of  structural  deformations  on  aeroelastic  loads.  This  situ¬ 
ation  required  detailed  aerodynamic  and  structural  analyses  which  were  computa¬ 
tionally  very  large.  In  some  of  the  studies  presented  in  reference  1,  composite 
materials  were  assumed  to  be  used  for  wing  construction.  Composite  materials 
allow  the  designer  the  added  flexibility  of  adjusting  the  proportions  of  mate¬ 
rial  at  the  various  orientations  and  the  orientation  angles  themselves.  This 
added  design  flexibility  increases  the  number  of  design  variables  and  further 
adds  to  the  complexity  of  the  sizing  problem. 

As  indicated  in  reference  4,  most  available  procedures  for  optimizing 
composite  structures  are  for  isolated  components  of  large  structures.  (See 
refs.  5  to  7.)  In  addition,  there  is  a  continuing  effort  for  improvement  in 
the  efficiency  of  optimization  procedures  to  extend  their  applicability  to 
entire  structures.  (See  refs.  8  to  13.)  Because  of  computational  size  and 
need  to  optimize  at  the  detail  design  level,  the  subject  problem  entails  char¬ 
acteristics  of  the  synthesis  on  both  component  and  large  structure  levels. 

The  purpose  of  this  paper  is  to  document  the  procedure  used  in  the  studies 
of  reference  1,  to  explain  the  algorithm  details,  and  to  give  a  more  detailed 
account  of  several  structural  sizing  trade-offs. 


Nick  J.  Santoro,  Southern  Company  Services,  Inc.,  Birmingham,  Alabama, 
participated  in  the  conceptual  part  of  the  reported  work  while  with  Vought  Cor 
poration,  Hampton  Technical  Center,  Hampton,  Virginia,  and  was  responsible  for 
generation  of  the  numerical  results. 


SYMBOLS 

Values  are  given  in  both  SI  and  U.S.  Customary  Units.  The  calculations 
were  made  in  U.S.  Customary  Units. 

A,B,D  matrices  of  membrane  stiffness,  bending  plate  stiffness,  and  bending 
membrane  coupling  stiffness,  respectively 

a,b  rectangular  panel  dimensions  (see  fig.  2) 

9i  dimensionless  measure  of  constraint  satisfaction 

I  set  of  constraints 

k  vector  of  plate  curvatures  due  to  bending 

M  vector  of  moments 

m  mass  of  core  face-sheet  bonding  material  per  unit  area 

mdep'mindep  nonoptimum  structural  mass  dependent  on  panel  area  and  indepen¬ 
dent  of  panel  area,  respectively 

mp  structural  mass  of  panel  defined  by  equation  (10) 

N  vector  of  membrane  force 

NX/Ny,NXy  normal  and  shear  membrane  forces  acting  on  edges  of  panel  (see 

fig.  2) 

S  panel  area 

t  thickness  of  material  whose  fibers  are  oriented  at  a  given  angle 

tc  "caliper"  thickness  of  sandwich  panel  (face-to-face) 

tt  sum  of  thicknesses  of  both  face  sheets 

t^jty  thickness  of  material  whose  fibers  are  oriented  at  <j>  and  y , 
respectively 

ti'n  thickness  of  ith  ply  and  total  number  of  plies  in  one  face  sheet  of 

symmetric  sandwich 

y  behavior  variable 
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ya  limit  on  behavior  variable  (allowable  value) 

Z  defined  by  equation  (6) 

e  vector  of  membrane  strains  gx,  Gy,  GXy  or  g-j  ,  £ 2 r  e12 

p  mass  density  of  ply  material 

pc  mass  density  of  core 

a- 1  stress  in  fiber 

a 2  stress  in  matrix,  normal  to  cr-j  f  in  ply  plane 

O' |t  allowable  tension  limit  of  fiber 

a2t'a2c  allowable  tension  and  compression  matrix  limits 
T] 2  in-plane  shear  stress 

<j) orientation  angles  of  fibers  (see  fig.  2) 

Subscr iptss 

a  allowable  limit 

cr  critical 

Configurations : 

B  metal  baseline  construction 

LM/H  low  Young's  modulus/hybrid  construction 

HM/H  high  Young's  modulus/hybrid  construction 

LM/P  low  Young's  modulus/pure  composite  construction 

HM/P  high  Young's  modulus/pure  composite  construction 

LM/H-CA  conservative  allowable  strain  variant  of  LM/H  configuration 

OVERALL  DESIGN  APPROACH 

Airplane  flexibility  introduces  a  coupling  of  aerodynamic  loads  and  struc¬ 
tural  displacements.  This  coupling  can  be  accounted  for  by  organizing  the  anal¬ 
yses  involved  into  an  iterative  process.  The  flow  chart  in  figure  1  (from 
ref.  14)  depicts  a  particular  organization  of  the  iterative  process  used  in 
this  study. 
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Two  iterative  loops  are  shown  in  the  flow  chart  of  figure  1 .  Loop  I  begins 
with  the  computation  of  aerodynamic  loads  on  the  entire  airplane  whose  shape  is 
preset  to  the  defined  cruise  shape  (best  for  cruise  aerodynamic  efficiency) . 
Initially,  the  same  shape  is  used  to  compute  loads  for  both  cruise  and  maneuver. 
The  loads  are  subsequently  applied  to  the  airframe  using  the  initial  cross- 
sectional  dimensions,  and  the  ensuing  analysis  yields  displacement  and  stresses 
for  both  cruise  and  maneuver.  The  cruise  displacement  outputs  are  used  to  cal¬ 
culate  the  jig  shape  (shape  to  which  the  airframe  is  to  be  built  to  attain  the 
desired  shape  under  cruise  loads).  The  jig  shape  and  the  maneuver  deflections 
are  superimposed  to  calculate  maneuver  shape.  The  loop  is  closed  by  feeding  back 
the  maneuver  shape  to  calculate  a  new  set  of  aerodynamic  loads  for  maneuver. 

In  the  flow  chart  of  figure  1 ,  loop  II  involves  use  of  the  stress  outputs 
to  resize  the  airframe  structural  components.  Resizing  of  structural  components 
produces  new  cross-sectional  dimensions  (e.g.,  wing-cover  thicknesses)  which  are 
fed  back  to  the  structural  analysis  input.  Following  the  method  established  in 
reference  14,  the  iteration  process  alternates  between  loops  I  and  II.  The  flex¬ 
ibility  of  the  system  allows  the  user  to  choose  the  execution  sequencing  that 
best  fits  the  problem  at  hand.  A  typical  mix  of  the  loops  is  three  executions 
of  loop  I  for  each  execution  of  loop  II.  Details  of  the  system  organization  are 
described  in  reference  15,  and  an  outline  of  the  major  analysis  programs  involved 
is  provided  in  reference  16.  All  data  transfers  showed  by  the  flow  chart  in 
figure  1  are  fully  computerized.  The  iterative  process  continues  until  the  con¬ 
vergence  criteria  for  aerodynamic  loads,  structural  mass,  and  stresses  are 
satisfied. 

Structural  resizing  is  executed  for  each  structural  component  separately 
as  if  that  component,  together  with  the  forces  acting  on  it  from  the  neighbor¬ 
ing  components,  were  extracted  from  the  surrounding  structure.  In  this  way,  an 
independent  problem  of  structural  optimization  is  formulated  with  respect  to  the 
component  forces  which  are  treated  as  invariant  external  loads.  Mathematical 
nonlinear  programing  is  used  to  solve  this  problem  for  each  structural  component. 
After  all  components  have  been  optimized,  new  aeroelastic  loads  are  computed, 
the  total  structure  is  reanalyzed,  and  new  component  forces  are  generated  for 
use  in  the  next  round  of  individual  component  optimizations. 

As  pointed  out  in  reference  17,  where  this  approach  and  its  positive  and 
negative  implications  are  discussed  in  detail,  the  disadvantage  of  such  a  decom¬ 
position  of  the  problem  is  that  optimization  of  components  is  not  tantamount  to 
optimization  of  the  whole;  or,  in  other  words,  minimization  of  the  individual 
component  masses  does  not  guarantee  minimization  of  the  total  mass.  This  situ¬ 
ation  is  caused  by  the  inability  to  control  the  load  path  on  the  assembled  struc¬ 
ture  level?  in  this  regard,  the  procedure  resembles  the  fully  stressed  design 
(FSD)  algorithm.  The  advantages  lie  in  the  drastically  reduced  dimensionality 
of  each  of  the  individual  optimization  problems  that  must  be  solved  and  in  the 
inherent  modularity  of  the  solution  algorithm.  The  consequence  of  the  former 
is  that  the  required  computer  resources  stay  within  reasonable  bounds.  The  mod¬ 
ularity  permits  selection  of  the  resizing  algorithm  considering  generality,  effi¬ 
ciency,  and  accuracy  in  isolation  from  the  rest  of  the  problem.  Details  of  the 
component  optimization  procedure  are  described  in  the  following  section. 
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OPTIMIZATION  PROCEDURE  FOR  SANDWICH  PANELS 

The  object  of  the  optimization  procedure  for  each  individual  panel  is  to 
minimize  the  panel  mass  subject  to  constraints  of  stress,  strain,  panel  buckling, 
and  minimum  gage.  Panel  load  consists  of  the  membrane  forces  Nx,  Ny,  and  Nxy 
shown  in  figure  2.  These  forces  are  held  constant  for  the  duration  of  each  con 
secutive  optimization  procedure  for  individual  panels.  For  metal  panels,  the 
design  variables  are  the  face-sheet  thickness  and  the  total  sandwich  thickness. 

For  composite  panels,  the  variables  are  the  thickness  of  the  composite  material 
laid  at  a  given  orientation  angle  (layer  thickness) ,  the  orientation  angles  them¬ 
selves,  and  the  total  sandwich  thickness.  Construction  of  a  composite  sandwich 
panel  and  the  design  variables  are  shown  in  figure  2.  Organization  of  the  pro¬ 
cedure  is  depicted  in  figure  3.  The  flow  chart  of  figure  3  corresponds  to  the 
box  marked  "RESIZE  ELEMENTS"  in  figure  1 .  A  panel  defined  by  the  design  vari¬ 
ables  and  subject  to  the  given  loads  is  analyzed  (box  1  in  fig.  3)  for  strain, 
stress,  and  panel  buckling.  Output  from  the  analysis  is  used  to  evaluate  con¬ 
straints  and  the  objective  function  which  is  the  mass  of  the  panel.  Subsequently, 
the  results  of  this  evaluation  are  used  by  a  general  purpose  optimizer  (box  2) 
to  calculate  a  new  set  of  design  variables  modified  to  reduce  the  objective  func¬ 
tion  and  to  move  toward  satisfaction  of  the  panel  constraints.  Operations  1 
and  2  are  repeated  until  the  constrained  minimum  mass  is  reached.  Then,  the  pro¬ 
cedure  is  repeated  for  each  next  consecutive  panel  until  all  panels  have  been 
optimized.  When  all  panels  have  new  dimensions,  the  wing  structure  is  reanalyzed 
for  aerodynamic  loads  and  internal  forces  (loops  I  and  II  in  fig.  1)  and  another 
round  of  the  individual  panel  optimization  procedure  is  carried  out. 


Evaluation  of  Constraints 

The  optimization  algorithm  requires  that  the  degree  of  violation  or  satis¬ 
faction  of  the  constraints  be  stated  in  the  following  standard,  dimensionless 
form: 


9i(y/Ya)  =  0 


(i  e  I)  (1) 


where  g^  =  —  -  1  symbolizes  the  ith  constraint,  y  is  a  behavior  variable 
Ya 

such  as  stress,  etc.,  and  ya  is  the  corresponding  allowable  value  of  y  dic¬ 
tated  by  the  material  limits  (yield  stress,  ultimate  stress,  etc.)  or  resulting 
from  the  structural  stability  considerations,  i.e.,  computation  of  the  critical 
buckling  stresses.  The  purpose  of  the  analysis  of  the  panel  being  optimized  is 
to  supply  values  of  y  and  ya  used  in  equation  (1). 

Side  constraints.-  The  side  constraint  (minimum  gage)  values  for  the  skins 
are  set  to  provide,  as  a  minimum,  four  plies  in  each  face  sheet  of  the  sandwich. 
These  plies  are  associated  in  pairs  with  orientation  angles  <t>  and  -<t>.  In 
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most  cases  <t>  is  selected  to  be  45°.  Thicknesses  associated  with  other  orien¬ 
tations,  that  is  0°  and  90°,  are  permitted  to  vanish  but  not  to  become  negative; 
therefore,  their  minimum  gage  is  zero. 

Upper  limit  on  the  skin  thicknesses  is  provided  indirectly  by  imposing  a 
minimum  gage  value  on  the  sandwich  core  thickness  as  a  fraction  of  total  sandwich 
thickness.  Without  such  constraint,  the  face  thicknesses  could  grow  inward  to 
the  extent  of  eliminating  the  core. 

Stress  constraints.-  Standard,  anisotropic  membrane  analysis  (ref.  18)  is 
used  to  produce  the  stress-strain  output  for  each  ply  for  given  edge  forces  and 
plate  geometry.  The  analysis  is  simplified  by  the  following  assumptions: 

(1)  The  face-sheet  plies  are  symmetric  with  respect  to  the  midsurface  of 
the  sandwich  panel. 

(2)  The  material  associated  with  each  orientation  is  uniformly  dispersed 
throughout  each  face  of  the  sandwich. 

(3)  The  lay-up  is  orthotropic  (t^  =  t^)  although  the  axes  of  orthotropy 
need  not  be  aligned  with  the  panel  sides  (general  orthotropy) .  As  shown  in 
figure  2,  the  orthotropy  axes  are  oriented  at  an  angle  y  with  respect  to  the 
x-axis  coinciding  with  the  panel  side. 

The  result  of  assumption  1  is  a  decoupling  of  the  basic  force-deformation  equa¬ 
tion,  from  reference  18, 


N 

A  B 

e 

M 

B  D 

k 

into  the  equations  N  =  Ae  and  M  *  Dk.  Assumption  2  permits  evaluation  of 
the  terms  of  matrices  A  and  D  independent  of  the  stacking  sequence,  and 
assumption  3  simplifies  the  buckling  analysis. 

Solution  of  the  membrane  problem  represented  by  N  =  Ae  yields,  for  each 
fiber  orientation,  stresses  parallel  and  perpendicular  to  the  fibers  a-j  and 
o 2r  the  in-plane  shear  stress  T 12 r  and  the  corresponding  strains  e*j  ,  e 2 / 
and  £]2-  These  values  are  combined  with  the  allowable  limits  in  constraint 
equations  in  the  standard  form  shown  in  equation  (1).  For  example,  if  a-j  ^  0, 


al 

- 1=0 

ctait 


(3) 
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where  O' ^  is  an  allowable  tension  and  C^.  is  an  associated  factor  used  to 
modify  the  allowable  stress  for  selected  load  cases.  A  typical  example  is  the 
reduction  of  the  allowable  stress  level  for  a  lg  cruise  condition  due  to  fatigue 
considerations.  All  numerical  information  on  load  factors  and  allowable  stress 
reduction  factors  actually  used  in  this  study  is  reported  in  the  section  on  the 
vehicle  model. 


In  order  to  prevent  matrix  failure,  the  strength  constaints  are  of  the  form 
given  in  reference  19: 


(4) 


where  a2,  =  a2t  if  a2  £  0,  or  a2a  =  02c  if  o2  <  0.  This  constraint  (eg.  (4)) 
is  a  dimensionless  inequality  equivalent  to  the  standard  form  represented  by  equa¬ 
tion  (1).  All  stress  constraints  are  evaluated  for  each  fiber  orientation  and  for 
each  loading  case.  In  addition  to  these  constraints,  strain  components  of  a  sand¬ 
wich  panel  are  kept  from  exceeding  specified  limits  by  setting  the  following  con¬ 
straints: 


x,a 


where  the  subscript  x  (or  y)  refers  to  the  panel  coordinate  that  has  at  least 
one  axis  aligned  with  one  of  the  panel  sides  and  the  subscript  a  denotes  corre¬ 
sponding  allowable  strains.  These  overall  constraints  are  useful  since  they  pro¬ 
vide  a  degree  of  additional  control  over  the  panel  strain  levels.  This  control 
is  important  for  at  least  two  reasons:  (1)  To  meet  the  structure's  overall  stiff¬ 
ness  requirements,  i.e.,  those  dictated  by  anticipated  aeroelastic  problems  while 
still  at  the  strength  design  stage  and  (2)  To  protect  the  panel  edge  members  from 
over stress  which  may  occur,  for  example,  in  the  case  of  titanium  spar  caps  fram¬ 
ing  a  composite  panel.  (See  section  "Types  of  Construction.”) 

Bucklinq  constraints.—  Buckling  of  'a  plate  such  as  the  one  shown  in  figure  2 
cannot  be  solved  exactly  in  a  closed  analytical  form  because  of  the  general  ortho- 
tropy  and  arbitrary  in-plane  loads.  Since  numerical  buckling  analysis  is  too 
costly  in  the  repetitive  application  required  in  optimization,  an  approximate 
but  fast  method  is  used  involving  analytical  solutions  given  in  reference  20  for 
plates  of  special  orthotropy.  The  special  orthotropy  assumption  results  in 
neglecting  the  bending- twist  coupling  terms  D16  and  D26  and  appears  to  be 
justified  for  this  type  of  application  by  the  findings  of  reference  21 . 
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According  to  reference  20,  the  general  state  of  load  is  separated  into  the 
cases  of  biaxial  normal  stress  and  shear  stress.  The  measure  of  buckling  criti¬ 
cality  £ue  to  biaxial  membrane  forces  is  given  by  the  quantity 


where  m  and  n  are  the  number  of  half-waves  of  the  buckling  deformation  func¬ 
tions  in  the  x-  and  y-directions ,  respectively. 

The  plate  buckles  when  Z  M ,  for  any  combination  of  m  and  n.  The  com¬ 
bination  that  minimizes  Z  is  searched  for  numerically.  A  conservative  approx¬ 
imation  to  the  critical  shear  is  obtained  by  the  formula  (ref.  20) ,  for  an  infi¬ 
nite  length  strip, 

[D1 1  <D22)  3]1//4 

(Nxy)cr  =  4c  “  (7) 

b2 
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Objective  Function 


The  structural  mass  of  the  panel  mp  is  the  objective  function.  The  vari¬ 
able  part  of  that  mass  is  the  sum  of  face-sheet  mass  and  core  mass  and  is  a  func 
tion  of  the  thickness  variables.  In  addition,  a  constant  part  is  also  included 
to  represent  masses  of  the  bonding  and  of  the  nonstructural  material  (m^p  and 
mindep)  according  to  the  formula 

n 

mp  =  2^~|j)tj_S  +  pcS(tc  -  tt)  +  mS  +  m^pS  +  n>indep[j  (10) 

i=1 


Optimization  Algorithm 


The  optimization  algorithm,  represented  by  box  2  in  figure  3,  used  for  the 
optimization  of  composite  panels  is  a  method  of  feasible-usable  directions,  orig¬ 
inally  due  to  Zoutendi jk  (ref.  22),  coded  in  a  program  described  in  reference  23. 
This  algorithm  turned  out  to  be  efficient,  converging  relatively  independent  o 
the  number  of  design  variables  for  the  reported  application. 


VEHICLE  MODEL 

The  structural  sizing  algorithm  described  in  the  foregoing  section  is  applied 
to  the  structural  box  of  the  low-aspect-ratio  arrow  wing  of  a  super sonic- transport 
configuration.  The  basic  details  of  the  finite-element  model  of  that  configuration, 
as  well  as  loads,  materials,  and  design  variables,  are  given  in  the  following 
section. 


Finite-Element  Model 

The  finite-element  model  is  shown  in  figure  4.  It  is  symmetric  with  respect 
to  the  airplane  longitudinal  center  line  with  one-half  consisting  of  746  grid 
points  and  2141  degrees  of  freedom  connected  by  2369  structural  finite  elements. 
The  model  represents  wing,  fuselage,  vertical  and  horizontal  tails,  wing  fins, 
engines,  and  engine  mounts.  Nonstructural  masses  of  the  engines,  leading-  and 
trailing-edge  devices,  and  fuel  tanks  with  fuel  are  included.  The  wing  is 
modeled  as  follows:  . 

(1)  Spar  and  rib  caps  -  rod  elements  having  only  ex tens ional  stiffness 

(2)  Spar  and  rib  webs  -  quadrilateral  panels  having  only  shear  stiffness 

(3)  Wing  covers  -  quadrilateral  and  triangular  membrane  elements,  isotropic 
for  metal  covers  and  anisotropic  for  composite  covers 

(4)  Engine  mounts  -  bar  elements  having  extensional,  bending,  and  torsional 
stiffnesses 
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(5)  Wing  fins  and  tail  -  quadrilateral  or  triangular  plane  elements  having 
in-plane  and  bending  stiffnesses 

(6)  Fuselage  frames  -  rod  elements 

(7)  Fuselage  longerons  -  rod  elements 

(8)  Fuselage  skin  -  quadrilateral  or  triangular  membrane  elements 

Since  the  fuselage  structure  is  not  included  in  the  resizing  process,  the  finite- 
element  model  of  the  fuselage  is  less  detailed  than  the  wing.  A  box  beam  is  used 
for  the  fuselage  and  is  composed  of  quadrilateral  membrane  elements  and  shear 
panels  preserving  the  mass  distribution  of  the  fuselage  and  its  overall  bending 
and  torsional  stiffnesses.  Statistically  derived  distributed  nonoptimum  and 
nonstructural  masses  are  accounted  for  by  appropriate  coefficients. 


Loading  Cases 

From  the  multitude  of  loading  cases  usually  considered  in  airframe  design, 
the  three  cases  selected  and  consistently  used  in  this  study  were:  2.5g  maneuver 
at  Mach  1.2,  -2g  taxi  at  Mach  2.7,  and  lg  cruise  at  Mach  2.7.  These  cases  consti¬ 
tute  a  minimal  set  of  loadings  needed  for  comparisons  of  the  construction  alterna¬ 
tives  and  for  methodology  evaluation.  The  rationale  for  their  selection  is  as 
follows:  (1)  the  maneuver  loads  generate  the  largest  bending  and  shear  forces 

near  the  wing  root;  (2)  negative  taxi  loads  activate  buckling  constraints  in  the 
bottom  cover  of  the  wing?  and  (3)  constraining  the  level  of  stress  due  to  the 
cruise  loads  is  a  simplified  way  to  control  airframe  fatigue  life.  Because  of 
aircraft  flexibility  and  difference  in  Mach  numbers,  the  maneuver  and  cruise  load 
distributions  differ  significantly  from  each  other.  Therefore,  stresses  corre¬ 
sponding  to  one  of  these  two  load  cases  cannot  be  obtained  by  scaling  stresses 
caused  by  the  other . 


Types  of  Construction 

In  this  study,  a  structure  made  entirely  of  titanium  is  adopted  as  a  base¬ 
line  design  against  which  the  composite  construction  results  are  compared.  Two 
types  of  composite  construction  are  considered  for  the  wing:  (!)  a  mixed,  or 
hybrid,  type  in  which  only  the  sandwich-panel  face  sheets  are  replaced  with 
composite  material  while  the  remainder  of  the  structure  is  made  of  metal  and 
(2)  a  pure  composite  construction  in  which  all  parts,  except  the  aluminum  honey¬ 
comb  sandwich  core,  are  made  of  composite  material. 


Material 

The  titanium  data  used  are  displayed  in  table  I.  The  allowable  stress 
level  for  cruise  is  restricted  to  the  value  of  Fac  to  approximately  account 
for  fatigue  requirements.  In  titanium  construction,  this  value  corresponds  to 
a  notch  factor  of  4.0.  The  data  for  the  honeycomb  core  and  core  face-sheet 
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bonding  are  also  included  in  table  I.  The  composite  material  used  in  hybrid 
and  pure  composite  construction  is  graphite-polyimide. 

Graphite-polyimide  is  a  composite  material  of  interest  for  supersonic- 
cruise  applications  because  of  its  relatively  good  retention  of  mechnical 
properties  at  elevated  temperatures.  The  material  properties  anticipated  for 
1986  (ref.  2)  are  displayed  in  table  I.  In  addition,  a  reduced  allowable  strain 
is  considered  to  establish  a  safe  bound  on  the  fiber  stress  allowables.  The 
fiber  volume  is  assumed  to  be  60  percent  throughout.  The  composite  is  assumed 
available  in  a  low  Young's  modulus,  high-strength  version  and  in  a  high  Young's 
modulus,  low-strength  version.  In  the  discussion  to  follow,  use  of  the  low- 
modulus  version  is  assumed  unl.ess  otherwise  indicated.  Insofar  as  the  fatigue 
stress  allowable  for  cruise  is  concerned,  no  data  similar  to  those  used  for 
titanium  construction  exist  for  composites.  However,  for  consistency,  the  same 
ratio  of  the  allowable  cruise  stress  to  the  design  limit  stress  is  used  for  com¬ 
posites  and  for  titanium. 


Design  Variables  and  Fixed  Parameters 

In  the  design  studies  reported  in  the  next  section,  wing  area,  planform, 
airplane  mass,  and  structural  layout  are  kept  constant.  Subject  to  change  are 
the  type  of  construction  and  the  material.  Within  each  of  these  discrete  choices, 
the  continuous  variables  for  hybrid  and  pure  composite  constructions  are  selected 
from  the  set  and  include  sandwich  depth  ("caliper"  thickness)  tc;  thicknesses  of 
the  layers  tg/  U ,  and  tgg  associated  with  orientation  angles  0°,  ±4>»  and 
90°,  respectively;  and  the  orientation  angles  y  and  <}> .  (See  fig.  2.)  Thus, 
the  largest  number  of  variables  is  six  since  the  balanced  laminate  assumption 
of  t.  =  t_A  holds.  In  order  to  address  the  obvious  manufacturing  requirements, 
grouping  of  some  of  these  design  variables  over  blocks  of  adjacent  panels  is 
also  considered.  For  the  case  of  metal  construction  (baseline) ,  the  variables 
are  sandwich  thickness  tc  and  face  thickness  tt. 


RESULTS 

Numerical  results  obtained  for  the  arrow-wing  structure  are  presented  in 
this  section.  The  results  pertain  to  three  major  consecutive  phases  of  the 
study:  (1)  comparison  of  construction  type  and  material,  (2)  optimization  with 

the  design  variables  determined  most  effective,  (3)  reoptimization  with  reduced 
number  of  design  variables  (grouping)  to  assess  penalty  for  uniformity  required 
for  manufacturing.  Results  discussed  in  this  report  represent  a  strength  design 
in  which  no  attempt  is  made  to  control  the  wing  stiffness  and  flutter. 


Comparison  of  Construction  Type  and  Material 

Several  alternative  designs  are  compared  in  terms  of  total  sandwich-cover 
structural  masses,  distributions  of  the  sandwich-cover  face-sheet  thickness  over 
the  wing,  criticality  of  the  buckling  constraints,  and  stiffness  characteristics 
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Thickness  distribution.-  The  distribution  of  total  thickness  tt  is  shown 
for  the  titanium  baseline  configuration  in  figure  5(a)  and  for  the  composite 
configurations  in  figures  5(b),  (c) ,  (d) ,  and  (e)  .  The  corresponding  masses 
are  recorded  in  table  II.  The  masses  represent  the  total  of  two  wings,  includ¬ 
ing  nonoptimum  mass. 

In  this  phase  of  the  study,  the  design  variables  are  face-sheet  thickness 
for  the  titanium  baseline  configuration  and  thicknesses  associated  with  0°,  ±45°, 
and  90°  angles  of  fiber  orientation  for  the  composite  configurations,  with  sand¬ 
wich  depth  fixed  at  tc  =  1  in.  for  all  configurations.  A  high  Young's  modulus, 
low-strength  version  of  the  composite  material  is  used.  As  seen  in  figure  5(a), 
the  distribution  is  characterized  by  the  thickness  plateaus  A  and  B,  which  corre¬ 
spond  to  the  location  of  the  wing  fin  and  engine  mount,  respectively.  A  similar 
thickness  distribution  for  the  hybrid  construction  can  be  seen  in  figure  5(b). 
Maximum  thicknesses  for  the  hybrid  construction  (figs.  5(b)  and  (d) )  are  about 
twice  the  thickness  for  titanium  construction.  However,  this  is  more  than  com¬ 
pensated  for  by  the  lower  density  of  the  composite  material,  which  results  in 
the  mass  savings  discussed  subsequently. 

Eliminating  the  titanium  caps  and  changing  the  shear  web  material  from 
titanium  to  ±45°  composite  turns  the  structure  into  a  pure  composite  construc¬ 
tion  with  the  skin  thickness  distribution  shown  in  figure  5(c).  With  titanium 
caps  removed,  the  composite  must  carry  a  greater  load;  therefore,  thickness  builds 
up  (as  seen  in  fig.  5(c))  to  values  higher  by  about  25  to  50  percent,  compared 
with  the  hybrid  construction. 

Thickness  contours  of  the  individual  plies  of  90°  and  ±45°  orientation  angles 
are  shown  in  figure  5(e)  for  the  LM/H  configuration.  The  contours  (fig.  5(e)) 
indicate  a  characteristic  concentration  of  the  ±45°  orientation  material  in  the 
area  where  the  load  path  curves  toward  the  wing  root  from  the  direction  established 
by  the  wing  sweep  angle.  Concentrations  typical  for  total  skin  thickness  are 
reflected  in  the  distribution  of  the  90°  orientation  (spanwise)  material. 

In  some  instances  these  concentrations  illustrate  a  well-known  instability 
potentially  inherent  in  a  fully  stressed  design  approach.  Under  that  approach, 
an  overstressed  element  is  stiffened,  thus  attracting  more  load  which  in  turn 
requires  more  stiffness  in  the  next  resizing  cycle.  An  example  of  such  thickness 
buildup  is  designated  "region  C"  in  figure  5(d).  This  numerical  instability  can 
be  checked  by  imposing  upper-limit  side  constraints  on  the  thickness  variables 
in  the  optimizer. 

Criticality  of  buckling  constraints.-  The  panels  which  evolve  from  the 
optimization  as  being  critical  in  buckling  are  shown  in  figures  6(a)  and  (b)  . 
Comparison  with  figure  5(d)  indicates  correspondence  of  the  thickness  ridge  and 
row  A  of  buckling  critical  panels  in  figure  6(a).  Removal  of  the  titanium  spar 
caps  puts  all  the  load  on  the  composite  panels,  thus  increasing  very  signifi¬ 
cantly  the  buckling  critical  area  of  a  low  Young's  modulus  wing  cover.  (See 
fig.  6(b)  .)  The  change  to  high  Young's  modulus  material  is  effective  in  reduc¬ 
ing  the  number  of  buckling  critical  panels  to  just  two  panels  in  figure  6(a) 
and  one  panel  in  figure  6 (b) . 
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Structural  mass  and  stiffness  characteristics.-  Attention  is  now  turned 
to  such  overall  structural  characteristics  as  structural  mass  and  stiffness. 

Total  structural  mass  of  two  wings  includes  mass  of  face  sheets,  core,  bonding, 
ribs,  and  spars.  Both  webs  and  caps  are  accounted  for  in  the  latter,  except 
for  a  pure  composite  construction  in  which  there  are  no  structural  caps.  Appro¬ 
priate  nonoptimum  mass  is  added  to  all  structural  components.  Static  stiffness 
can  be  characterized  by  wing  tip  deflection  (relative  to  the  wing  root)  under 
2.5g  pull-up  maneuver  load,  and  stiffness  can  be  partially  measured  by  the 
natural  frequency  of  the  wing  first  bending  mode. 

The  three  wing  characteristics,  mass,  tip  deflection,  and  fundamental  fre¬ 
quency,  are  represented  by  a  vertical  bar  diagram  in  figure  7  for  the  B,  LM/H, 
HM/H,  LM/P,  HM/P,  and  LM/H-CA  configurations.  Comparison  of  the  bars  shows  a 
significant,  20-percent  decrease  in  structural  mass  brought  about  by  changing 
from  the  B  configuration  to  the  LM/H  configuration.  That  benefit  is  partly  off¬ 
set  by  an  increase  in  the  wing  tip  deflection  and  a  loss  in  the  frequency  of  the 
wing  first  bending  mode. 

The  loss  of  stiffness  is  caused  by  the  fact  that  the  ratio  of  the  allowable 
stress  of  the  composite  to  that  of  titanium  is  higher  than  the  ratio  of  the 
respective  Young's  moduli,  and  is  also  strongly  influenced  by  changes  in  thick¬ 
nesses  associated  with  different  orientation  angles.  Lower  strength  but  higher 
Young's  modulus  of  the  HM/H  configuration  reduces  the  mass  by  13  percent  and 
increases  stiffness,  as  manifested  by  a  drop  in  the  wing  tip  deflection. 

A  drastic  mass  reduction  of  56  percent  is  observed  in  the  LM/P  configura¬ 
tion,  but  it  is  countered  by  an  equally  sharp  increase  of  the  wing  tip  deflection. 
However,  there  is  a  small  decrease  of  frequency  because  of  the  attendant  decrease 
in  mass.  These  effects  are  caused  by  the  removal  of  the  titanium  caps  with  their 
contribution  to  the  mass  and  their  limiting  influence  on  the  composite  strains. 
(See  previous  discussion  on  the  constraints  protecting  the  titanium  caps  from 
over stress  in  a  hybrid  construction.) 

The  HM/P  configuration  diminishes  the  mass  savings  to  49  percent  but  radi¬ 
cally  improves  stiffness,  as  indicated  by  significant  reduction  of  the  wing  tip 
deflection  and  rise  in  the  wing  frequency,  the  latter  being  tempered  somewhat 
by  growth  of  the  structural  mass.  The  influence  of  the  imposition  of  a  conserva¬ 
tive  strain  allowable  in  the  LM/H  variant  (m/H-CA  configuration)  is  seen  to  cause 
the  mass  reduction  to  drop  to  16  percent  (from  21  percent  noted  previously). 

The  wing  tip  deflection  decreases  slightly  from  the  baseline,  but  the  wing  fre¬ 
quency  increases  because  of  reduction  in  mass. 

In  summary,  data  collected  in  figure  7  and  table  II  show  the  maximum  mass 
savings  to  be  56  percent  of  the  baseline  mass  when  the  configuration  used  is  the 
low  Young's  modulus  pure  composite  construction  (LM/P). 


Selection  of  Design  Variables 

The  subset  of  design  variables  consisting  of  three  thicknesses  for  orienta¬ 
tion  angles  fixed  at  0°,  ±45°,  and  90°  has  been  used  so  far.  Effects  of  freeing 
the  other  design  variables  are  now  assessed  for  the  LM/H-CA  configuration. 
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Distributions  of  thicknesses  and  orientation  angles.-  Distributions  of  the 
thicknesses  are  illustrated  in  figures  8(a)  and  (b)  ,  which  show  contours  of  total 
thickness  of  all  layers  in  both  faces  of  the  upper  wing  cover  sandwich  panels. 
Figure  8(a)  corresponds  to  the  layer  thicknesses  tgf  t(j)f  tgg,  sandwich  depth 
tc,  and  orientation  angle  (j) .  All  thicknesses  t  and  tc  are  free  design  vari¬ 
ables  while  orientation  angle  y  is  kept  constant  at  a  value  of  90°.  Figure  8(b) 
corresponds  to  the  set  of  design  variables  reduced  to  tg,  ,  tgg,  and  y 

while  freezing  tc  and  (j)  (<J>  =  45°);  figure  8(c)  illustrates  the  result  of  free¬ 

ing  the  variable  tc  and  thicknesses  tg,  ,  and  tgg  while  freezing  the 
angles  y  and  (j)  at  90°  and  45°,  respectively;  and  finally,  figure  8(d)  depicts 
the  result  for  all  six  design  variables  being  simultaneously  free.  Comparison 
of  the  results  in  figure  8  shows  significant  changes  of  the  cumulative  ply  thick¬ 
ness  distribution  caused  primarily  by  the  smoothing  effect  of  the  variable  tc. 
This  effect  is  due  to  the  fact  that  the  sandwich  depth  increase  is  a  very  effi¬ 
cient  means  of  alleviating  the  panel  buckling. 

Individual  ply  thicknesses  and  orientation  angles  $  are  presented  in  fig¬ 
ure  9(a)  for  variable  $  and  constant  y.  Orientation  angle  is  indicated  in 
figure  9(a)  by  arrows  for  each  panel  (numerals  at  each  arrow  denote  the  associ¬ 
ated  ply  thickness) .  The  orientation-angle  distribution  reveals  that  diagonal 
fibers  associated  with  angle  <f>  tend  to  align  with  the  spanwise  fibers  in  the 
rear  of  the  inboard  wing  box,  which  is  dominated  by  bending  moment.  A  signifi¬ 
cant  number  of  angles  $  close  to  conventional  ±45°  orientation  is  found  in 
the  front  upper  portion  of  the  inboard  wing  box,  where  the  load  path  is  turning 
toward  the  fuselage  as  it  comes  from  the  outboard  wing. 

Numerical  values  of  angles  y  and  4>  when  both  are  used  as  optimization 
variables  are  shown  in  figure  9(b).  It  is  apparent  that  the  original  setting 
of  y  =  0°  (or  close  to  it)  is  preserved  by  the  optimizer  over  a  great  portion 
of  the  wing  area.  The  only  wing  areas  where  it  is  advantageous  to  alter  y 
in  the  optimization  process  are  in  the  vicinity  of  the  wheel  well,  over  the  out¬ 
board  wing  tip,  and  near  the  fin  attachment,  where  the  structural  layout  causes 
abrupt  changes  of  the  load  path  direction. 

Influence  on  buckling  behavior.-  The  location  of  the  buckling  critical 
panels  over  the  wing  depends  significantly  on  the  choice  of  design  variables, 
as  revealed  by  comparing  the  different  combinations  in  figures  10(a),  (b) ,  (c) , 
and  (d)  with  figure  6(a).  Adding  more  design  variables  permits  satisfaction 
of  the  stress-strain  constraints  with  less  material  and,  therefore,  makes  buck¬ 
ling  significantly  more  critical;  this  is  illustrated  by  the  greater  number  of 
buckling  critical  panels  in  figures  10(b)  and  (d)  than  in  figure  6(a).  Increase 
in  the  number  of  such  panels  over  the  lightly  loaded  front  part  of  the  inboard 
wing  (fig.  10(c))  indicates  that  the  sandwich  depth  tc  is  decreased  to  reduce 
the  core  mass. 

Changes  in  mass  and  stiffness  characteristics.-  Variations  of  the  thickness 
distributions  caused  by  different  choices  of  design  variables  have  a  cumulative 
impact  on  total  mass,  as  indicated  by  values  collected  in  table  III.  The  table 
is  arranged  to  show  progressively  lower  mass  from  top  to  bottom  brought  about  by 
the  design  variable  choices  listed.  This  arrangement  reveals  a  general  trend 
of  structural  mass  decrease  brought  about  by  greater  design  freedom,  that  is,  a 
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larger  number  of  design  variables  which  are  free  to  change  under  the  control  of 
the  optimization  algorithm.  Thus,  a  change  from  "thickness  only"  to  "all  six" 
design  variables  reduces  the  mass  of  the  sandwich  covers  (sum  of  column  3  and 
column  4)  by  21.5  percent  and  the  total  mass  by  9.4  percent.  A  comparison  of 
mass  differences  points  to  the  sandwich  depth  tc  as  the  most  effective  vari¬ 
able  in  reducing  mass.  Angle  <t>  consistently  causes  more  change  than  angle  y, 
which  is  the  least  effective  variable.  The  effect  of  4>  is  somewhat  magnified 
when  coupled  with  tc,  as  shown  by  comparison  of  relative  mass  decrements  in  the 
third  and  seventh  rows  of  table  III. 

The  wing— tip- def lection  column  of  table  III  shows  that  activating  the  angle 
design  variables  reduces  stiffness  more  than  structural  mass.  The  variable  tc 
has  a  similar  effect,  although  it  appears  to  be  less  detrimental  to  stiffness 
than  the  variables  <j>  and  y.  The  loss  of  stiffness  which  accompanies  the  mass 
reduction  can  be  attributed  to  the  lack  of  the  stiffness  constraints  in  the 
optimization. 

Results  presented  in  this  section  identify  the  sandwich  depth  tc  and  orien¬ 
tation  angle  4>  as  the  design  variables,  besides  the  ply  thicknesses,  which 
appear  to  be  most  effective  under  the  particular  set  of  loads  and  constraints. 
However,  there  are  obvious  manufacturing  constraints  that  limit  the  design  vari¬ 
able  freedom.  It  lies  within  the  procedure  capability  to  address,  at  least  par¬ 
tially,  some  of  these  constraints,  as  described  in  the  following  section. 


Variables  Constant  Over  Large  Areas  of  Wing  Cover 

In  the  foregoing  discussion,  the  manufacturing  constraints  are  set  aside, 
temporarily,  in  order  to  permit  an  unrestricted  exploration  of  the  benefits  that 
may  be  brought  about  by  freeing  those  design  variables,  such  as  orientation  angles, 
which  under  a  conventional  approach  would  probably  remain  fixed. 

Among  the  design  variable  choices  considered,  the  proposition  of  varying 
orientation  angles  $  and  y  for  each  panel  separately  is  the  least  producible 
one.  However,  inspection  of  the  <t>  and  y  distributions  over  the  wing,  such 
as  shown  in  figure  9(b),  suggests  the  possibility  of  grouping. 

Examples  of  grouping  for  variables  y  and  <t>  are  illustrated  by  fig¬ 
ures  11(a)  and  11(b),  respectively.  Similarly,  grouping  of  sandwich  depth  tc 
is  depicted  in  figure  11(c).  The  effect  is  that  rather  large  wing  areas,  e.g., 
almost  the  whole  "glove"  region,  may  be  assigned  a  single  value  of  a  particular 
variable.  The  question  of  whether  or  not  the  groupings  shown  go  far  enough  to 
render  the  extra  manufacturing  cost  justifiable  can  only  be  answered  on  a  case- 
by-case  basis;  the  point  to  be  made  here  is  simply  that  such  a  grouping  can  be 
made. 


Examples  of  two  grouping  alternatives  are  illustrated  in  figure  12  for 
both  upper  and  lower  wing  covers.  Grouping  2  (fig.  12(b))  is  more  restrictive 
than  grouping  1  (fig.  12(a))  as  it  results  in  a  more  uniform  (more  producible) 
distribution  of  angle  <t>.  Once  a  grouping  is  defined,  the  optimization  is 
repeated  with  the  grouped  variables  input  as  constants  in  the  optimization  pro¬ 
cedure.  The  purpose  of  this  operation  is  to  adjust  the  remaining  free  variables 
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to  the  new  situation  (new  load  path,  new  constraints  becoming  critical)  which 
is  created  by  altering  the  previous  optimization  result  by  grouping.  Such 
adjustment  by  means  of  reoptimization  restores  a  major  part  of  the  mass  reduc¬ 
tion  gained  by  means  of  the  original  optimization  (ref.  24) .  The  results  of 
the  reoptimization  are  shown  in  table  IV  in  the  form  of  sandwich  cover  face- 
sheet  total  mass  for  the  two  groupings  portrayed  in  figures  12(a)  and  (b) .  The 
table  lists  a  9.3-percent  mass  reduction  gained  by  freeing  4>  in  addition  to 
thicknesses,  when  compared  with  a  reference  (0°,  ±45°,  90°)  design  in  which  only 
the  thicknesses  are  free.  That  mass  reduction  diminishes  to  7.8  percent  and 
6.2  percent  for  groupings  1  and  2,  respectively.  The  results  indicate  that  two- 
thirds  of  the  weight  reduction  realized  by  freeing  the  variable  <J>  can  still  be 
recovered  wnen  that  variable  is  grouped  and  frozen.  As  expected,  a  more  restric¬ 
tive  grouping  corresponds  to  a  lower  percentage  of  recovery. 


Summary  of  Trends 


Results  reviewed  in  the  preceding  section  appear  to  have  design  implications 
which  can  be  summarized  as  follows. 

Structural  mass  savings  of  56  percent  can  be  obtained  for  the  subject  airframe 
sized  for  strength  only  by  changing  from  baseline  titanium  to  pure  composite  con¬ 
struction  of  a  high-strength  (low  Young's  modulus)  material.  Mass  reductions  due 
to  hybrid  construction  and  use  of  composite  material  having  lower  strength  but 
higher  Young's  modulus  values  are  contained  in  that  range.  These  mass  reductions 
correspond  to  the  use  of  a  conventional  0°,  ±45°,  90°  lay-up  and  a  constant  depth 
for  the  sandwich  wing  cover  panels. 

Further  mass  reductions  on  the  order  of  5  percent  are  possible  if  the  number 
of  design  variables  is  increased  to  include  orientation  angles  and  sandwich  depth. 
Most  of  that  increased  mass  reduction  potential  is  realizable  if  one  addresses 
the  producibility  constraints  by  reoptimization  with  grouping  and/or  linking. 

These  structural  mass  trends  are  depicted  in  figure  13,  which  shows  (from  left 
to  right)  the  structural  mass  of  the  construction  variants  discussed  in  the  fore¬ 
going,  relative  to  the  titanium  baseline  taken  as  100  percent. 

Except  for  the  thickness  design  variables,  sandwich  depth  consistently 
emerges  as  the  variable  most  effectively  influencing  structural  mass,  suggesting 
a  design  with  stepped  or  taper ed-wing  sandwich  cover  depths.  The  high  Young's 
modulus  material  option  appears  to  be  the  most  effective  means  to  counter  the 
stiffness  decrease  associated  with  the  change  to  composite  materials. 


Convergence  and  Computational  Requirements 

The  method's  convergence  is  illustrated  in  figure  14  by  iteration  histories 
for  the  two  configurations  B  and  HM/H.  Loops  I  and  II  (see  fig.  1)  are  typically 
performed  in  a  ratio  varying  from  1/1  to  3/1,  that  is,  one  to  three  aerodynamic 
loads  updates  (loop  I)  per  one  structural  resizing  (loop  II) .  One  complete  itera¬ 
tion  is  defined  as  one  execution  of  resizing  loop  described  by  the  flow  chart  in 
figure  1  plus  the  associated  load  update  loop  executions.  Quantities  plotted 
(fig.  14)  are  structural  mass,  wing  root  bending  moment,  and  wing  tip  deflection. 
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The  curves  indicate  a  smooth  monotonic  convergence  for  both  configurations. 

All  three  quantities,  and  particularly  the  root  bending  moment,  stabilize  very 
fast  in  the  baseline  configuration.  Convergence  is  not  as  fast  for  the  bending 
moment  and  tip  deflection  in  the  composite  construction  because  of  the  composite 
wing's  greater  flexibility,  and  consequent  stronger  coupling  between  the  struc 
tural  deflections  and  aerodynamic  loads.  The  mathematical  optimization  of  eac 
panel  of  the  wing  cover  requires  from  8  to  1 5  objective  function  and  constraint 
evaluations  and  terminates  mostly  by  the  criterion  of  small  consecutive  changes 
of  the  objective  function.  In  summary,  the  whole  iterative  procedure  appears 
to  be  well  behaved  and  converges  to  a  reasonable  tolerance  within  five  itera- 
tions.  Approximately,  this  consumes  as  much  computer  resources  as  seven  repeti¬ 
tions  of  a  complete  aerodynamic  load  and  finite-element  structural  analyses  of 
the  subject  airframe. 

The  procedure  depicted  in  figure  1  is  composed  of  sequential  execution  of 
several  computer  programs.  The  major  programs  are  listed  in  table  V  with  their 
field  lengths  and  CPU  time  for  one  execution  for  a  CONTROL  DATA  CYBER  175  com¬ 
puter  system.  Computer  performance  indicates  acceptable  consumption  of  computer 
resources  and  generally  produces  a  wing  optimization  within  five  iterations,  CPU 
700 -| o  sec,  and  less  than  24  hours  of  turnaround  time. 


CONCLUDING  REMARKS 

A  computerized  procedure  to  generate  cross-sectional  dimensions  and  fiber 
orientations  for  composite  airframe  structures  has  been  developed.  A  methodology 
of  its  application  to  a  wing  structural  synthesis  has  been  established  involving 
the  selection  of  the  type  of  construction,  material,  and  design  variables.  The 
algorithm  unifies  computations  of  aeroelastic  loads,  stresses,  and  deflections 
and  optimal  structural  sizing  and  fiber  orientations  in  an  open-ended  system  of 
integrated  computer  programs.  Aerodynamic  loads  analysis,  finite  element  struc¬ 
tural  analysis,  and  a  mathematical-optimization  technique  are  the  main  components 
of  the  procedure.  Each  wing  cover  panel  is  optimized  separately  for  minimum  mass 
subject  to  stress,  strain,  local  buckling,  and  minimum  gage  constraints.  The 
design  variables  are  ply  thicknesses  and  fiber  orientation  angles.  Experience 
gained  in  application  to  a  low-aspect-ratio  wing  of  a  large  supersonic  transport 
airframe  leads  to  the  following  conclusions  regarding  the  procedure: 

1 .  A  unified  iterative  method  combining  structural  resizing  and  aeroelastic 
load  computation  converges  smoothly  and  monotonically  for  both  metal  and  composite 

wings. 

2.  Use  of  mathematical  programing  in  an  integrated  open-ended  and  modular 
system  of  computer  programs  provides  flexibility  needed  for  variety  of  design 
variables,  types  of  material  and  construction,  and  constraints;  including  some 
of  the  producibility  requirements  that  have  to  be  dealt  with  in  composite  struc¬ 
tures.  However,  since  individual  panels  are  optimized  separately,  the  resulting 
design  is  expected  to  be  near,  but  not  at,  the  minimum  total  mass  design. 

3.  Computer  performance  of  the  method  as  implemented  on  the  CONTROL  DATA 
CYBER  175  computer  indicates  acceptable  consumption  of  computer  resources  and 


17 


generally  produces  a  wing  optimization  within  five  iterations,  CPU  700 -j q  sec, 
and  less  than  24  hours  of  turnaround  time. 

In  the  particular  application,  structural  mass  reductions  in  excess  of 
50  percent  result  when  changing  from  baseline  titanium  to  pure  composite  con¬ 
struction.  Only  about  one-third  of  that  mass  saving  is  realized  when  the  change 
is  made  to  a  hybrid  construction  which  retains  the  titanium  network  of  ribs  and 
spars. 

Varying  the  sandwich  panel  depth,  in  addition  to  the  layer  thicknesses,  is 
found  to  be  a  significantly  more  effective  means  of  reducing  the  structural  mass 
than  varying  the  ply  orientation  angles.  Generally,  the  mass  reduction  is  asso¬ 
ciated  with  loss  of  stiffness.  That  loss  can  be  reduced  by  use  of  a  higher 
Young’s  modulus  material. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
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TABLE  I.-  MATERIAL  PROPERTIES 


TABLE  II.-  MASS3  AND  STIFFNESS  DATA  FOR  SEVERAL  CONFIGURATIONS 


composite  configurations  are  for  fixed  core  depth  and  ply  orientations  of  0°f 


TABLE  III.-  FREE  DESIGN  VARIABLES  FOR  LM/H  CONFIGURATION3 
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TABLE  IV.-  RESULTS  OF  GROUPING  DESIGN  VARIABLE 
LM/H  CONFIGURATION 


TABLE  V.-  COMPUTER  RESOURCES 


Program 


Field  length 

CPU  time,  sec 

(a) 

SPAR  -  a  finite  element 

165  100 

analysis  program 

Jig-shape  computation 

40  100 

Woodward-Car michael 

52  100 

aerodynamics 

Generation  of  pressure 

52  100 

coefficients  for  camber 

and  twist 

Airplane-balance 

52  100 

calculation 

Rib  and  spar  web  resizing 

50  100 

Skin  resizing 

120  100 

L _ _ 

aCONTRO L  DATA  CYBER  175  computer  system 
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Figure  1  Iterative  procedures  for  aeroelastic  loads  computation  (loop  I)  and 
wing  cover  resizing  (loop  II)  (from  ref.  14). 


Composite  sandwich  panel  and  desi 


Figure  3.-  Resizing-procedure  flow  chart. 


2 


projections  of  finite-element  model.  Inset  shows  construction  details. 


Skin  thickness  contours  in 


Region 


Concluded. 


Buckling  critical  for  low  Young's  modulus 
Buckling  critical  for  high  and  low  Young's  modulus 


Structural  mass 
Wing  tip  deflection 
Frequency  of  wing  first  bendin 


thickness  contours  in  0.0254  cm  (0.01  in.)  for  different  combinations  of  design 
variables  for  LM/H  configuration. 


Buckling  critical  panels  for  different  combinations  of  design  variables 

for  LM/H  configurations. 
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(a)  Orientation  angle  y.  (Values  are  given  in  degrees.) 


(b)  Orientation  angle  <j>.  (Values  are  given  in  degrees.) 

Figure  11.-  Grouping  of  panels  for  design  variables  y,  c}) ,  and  tc. 
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Mass5  percent 


Figure  13*-  Effects  on  structural  mass  of  changes  in  construction  material, 

design  variables,  and  constraints. 
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